Phase-Sensitive Probes of Nuclear Polarization in Spin-Blockaded Transport 
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Spin-blockaded quantum dots provide a unique setting for studying nuclear-spin dynamics in a 
nanoscale system. Despite recent experimental progress, observing phase-sensitive phenomena in 
nuclear spin dynamics remains challenging. Here we point out that such a possibility opens up in 
the regime where hyperfine exchange directly competes with a purely electronic spin-flip mechanism 
such as the spin-orbital interaction. Interference between the two spin-flip processes, resulting 
from long-lived coherence of the nuclear-spin bath, modulates the electron-spin-flip rate, making it 
sensitive to the transverse component of nuclear polarization. In a system repeatedly swept through 
a singlet-triplet avoided crossing, nuclear precession is manifested in oscillations and sign reversal of 
the nuclear-spin pumping rate as a function of the waiting time between sweeps. This constitutes a 
purely electrical method for the detection of coherent nuclear-spin dynamics. 



Due to their long coherence times, nuclear spins offer 
unique opportunities to study quantum dynamics. While 
conventional techniques using RF fields probe macro- 
scopic groups of spins, currently there is wide interest in 
the local dynamics of nanoscale groups of spins. In par- 
ticular, semiconducting quantum dots have emerged as a 
platform for investigating the intriguing quantum many- 
body dynamics of coupled electron and nuclear spins [l|- 
,8]. Many experimental achievements in this field rely 
on the phenomenon of spin blockade, in which electrons 
must flip their spins in order to pass through the sys- 
tem in accordance with the Pauli exclusion principle [9( . 
Such spin-blockaded transport is sensitive to the nuclear- 
polarization-dependent Overhauser field, which controls 
the energy splittings between electronic spin states. This 
sensitivity affords an exciting opportunity to electrically 
control and detect the states of nuclear spins [ToL Hlj- 

In this Rapid Communication, we propose a method 
which can be used to electrically probe coherent nuclear 
phenomena, such as Larmor precession or Rabi oscilla- 
tions, which are not easily probed by the Overhauser ef- 
fect alone. We identify a regime where transport is sen- 
sitive to all three vector-components of the local nuclear 
polarization, due to coherent interplay between hyperfine 
coupling and an electron-only spin-coupling such as the 
spin-orbit interaction [13, [HI]. As we show, systems with 
such sensitivity host a variety of new phenomena aris- 
ing from coherent nuclear dynamics. Recently observed 
oscillations of the nuclear pumping rate controlled by nu- 
clear Larmor precession [FJ] indicate that this unexplored 
regime in now within experimental reach. 

Here we investigate these phenomena in a two-electron 
double quantum dot in a uniform magnetic field. The 
two-electron singlet and m — ±1 triplet states are cou- 
pled via the difference of hyperfine fields due to transverse 
nuclear polarization in the two dots, AB nuCj i 1^, 11 , ll| . 
Without an additional coupling, the electron spin-flip 
rate depends only on |AB nuCj jJ, and not on the orien- 
tation of AB nuc j. in the XY plane [lij ]. 

The situation becomes considerably more interesting 
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FIG. 1: (Color online) Coherent interplay of hyperfine and 
spin-orbit mediated transitions, a) Interdot tunneling in the 
triplet state is suppressed by Pauli exclusion, but can be me- 
diated by the hyperfine and spin-orbit interactions which do 
not conserve electron spin, b) The singlet and triplet lev- 
els \S) and | T+ } exhibit an avoided crossing with split- 
ting \vg\ — \vso + «HFe l9 | when exchange energy compen- 
sates Zeeman energy, c) Behavior near the S — T+ cross- 
ing, controlled by \vg\, is sensitive to the relative phase 8 
between spin-orbit and hyperfine matrix elements, d) Phase- 
dependent transition probability, Eq.@, with vso = OAy/Wfi, 
and vhf = 0.6-/H/3. 



when singlet-triplet transitions can occur due to either 
the spin-orbit or the hyperfine interaction (see FigJIJ . In 
this case, we find that the probability of electron spin- 
flip depends on the angle 9 of the nuclear polarization 
measured in the XY plane relative to a fixed axis deter- 
mined by the spin-orbit interaction. Such a dependence 
makes electron transport sensitive to the phase of nuclear 
spin precession. Below, we analyze this phenomenon 
for a Landau-Zener-type process, occurring when the 
electronic system is swept through a singlet-triplet level 
crossing, as employed e.g. in Refs.[14l. Il6|| and depicted 
in FigHJ). To make contact with experiment, we focus 
in particular on resulting signatures in nuclear spin po- 
larization. Strikingly, we find robust oscillations in the 
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nuclear spin pumping rate which result from nuclear pre- 
cession and survive even after averaging over the distri- 
bution of random initial nuclear spin states. Note that 
other types of electron-only coupling such as Zeeman cou- 
pling to the nonuniform field of a micromagnet 17] can 
produce analogous effects. 

The origin of the spin-angle dependence can be un- 
derstood heuristically by analyzing the avoided crossing 
that opens near the degeneracy of the triplet and singlet 
levels | T + ) and IS*) in nonzero magnetic field, circled 
in FigOJ). Due to the spin-orbit interaction, tunneling 
between dots is accompanied by a rotation of electron 
spin [li| which introduces a nonzero spin-flip amplitude. 
In addition, the hyperfine interaction between electron 
and nuclear spins allows electronic transitions accom- 
panied by nuclear spin flips, described by the effective 
Hamiltonian H KF = A+\ S)(T+ \ + A~\ T+ }( S |, with 

A ± = J^i^i" 9ilf ■ The magnitude and sign of each cou- 
pling constant gi depends on the location of nucleus i 
(positive in dot I and negative in dot 2). 

Taking advantage of the fact that the number of nu- 
clei interacting with the electrons is very large, iV nuc sa 
O(10 6 ), we note that the commutator L4 + ,A~] is typ- 
ically smaller than A + A~ by a factor of order 
This allows us to treat A + as a nearly classical variable, 
A + m i>hf e" 10 . Here uhf is proportional to the magnitude 
of the transverse hyperfine difference field AB nuc- j_, and 
9 describes its orientation in the XY plane. 

Together, the spin-orbit and hyperfine pieces provide 
a matrix element between the singlet and triplet states: 



Vg = ( S \H\ T + ) = vso + VHFe 1 



(1) 



Evolution in the {|<S), subspace near the level 

crossing is described by the 2x2 Hamiltonian 
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A(t)=s T+ (t)-e s (t), (2) 



where £t + and e$ are the energies of the diabatic | S ) and 
| T + ) states. The level detuning A(t) can be controlled 
by electrostatic gates and/or magnetic field. 

We consider the case where the system is initialized to 
the "(0,2)" singlet state \S) with both electrons resid- 
ing in the right dot (large positive A), and then swept 
through the avoided crossing to the "(1,1)" charge regime 
with one electron on each dot (see FigQ}j). For a con- 
stant sweep rate j3 = \dA/dt\, the electron spin flip in 
such a model is interpreted as a Landau-Zener transition 
occurring with probability (see FigQJl) 



Plz = 1 - cxp (-27T 
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where Plz is the probability to remain in the ground 
state. The explicit dependence on the phase 9 of nuclear 
polarization, which enters through the matrix element 
vg, shows the singlet /triplet transition probability's sen- 
sitivity to the transverse nuclear polarization vector. 



Model ([2]) provides a useful heuristic for understanding 
electron spin dynamics, in particular for situations where 
the nuclear spin state is characterized by a well-defined 
azimuthal angle 9. However, to understand the behav- 
ior with more general initial states and to account for 
the effects of back-action on the nuclei, we must examine 
the quantum many-body dynamics of coherently coupled 
electron and nuclear spins. By solving this problem be- 
low, we will further justify the form of Eq.([T]) and will 
obtain the electron and nuclear spin-flip rates. 

A key new feature here is the relaxed selection rule gov- 
erning nuclear pumping: due to spin-orbit coupling, elec- 
tron spin flips may occur with or without a compensating 
nuclear spin flip. Although these two processes appar- 
ently lead to orthogonal final states, long-lived coherence 
of the nuclear spin bath allows interference between dif- 
ferent transition sequences (see Figj2]). Moreover, a single 
electron can change the total nuclear polarization by an 
amount Am that can be larger than 1, and of either sign. 

At this point it is convenient to map the problem 
onto the bipartite 1-dimensional quantum walk shown 
in FigHK where each unit cell is labeled by to, the z- 
projection of the total nuclear spin I z = ^ If. In doing 
so, [J z ,j4 ± ] remains nonzero. Intracell hopping between 
internal states T and S, characterized by Am = 0, de- 
scribes a spin-orbit transition occurring with amplitude 
vso- Intercell hopping, characterized by Am = ±1, de- 
scribes a hyperfine transition that occurs with the ampli- 
tude «hf — (T, to— 1 | Phf \S, to). Initially, we consider 
sweeps which are fast compared to the nuclear Larmor 
period, and so neglect the nuclear Zeeman energy. 

Because N nuc is large, the value of uhf will change very 
little during the course of a few sweeps, and we treat it 
here as a constant. However, the initial values of whf 
may differ from one run to another, with a statistical 
distribution of the form p(v) oc ve~ v / s , where s is a 
constant. Additionally, because the typical values of total 
nuclear spin will be large, of order \fN nnc , we may take 
the ladder of allowed values of to to be infinite while the 
total number of sweeps is not too large. 

The state of the system evolves according to 
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At time t , the system is initialized to the singlet state 
with polarization m — m : ip^ = <5 m . mo , V'm = 0- The 
expected change of polarization, 
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is determined by the probabilities {P m } to find the nu- 
clear spin with z-projection to in the final state at time 
tp after the sweep. For convenience we take the (non- 
classical) initial nuclear state to be an eigenstate of J 2 , 
but the framework can be applied for more general states. 
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FIG. 2: (Color online) Quantum walk model of electron- 
nuclear spin dynamics near the S — T+ crossing, Eq.Q. a) 
Unit cells are labeled by m, the z-component of total nu- 
clear spin. Intracell and intercell hopping correspond to spin- 
orbit and hyperfine transitions, with matrix elements vso and 
dhf, respectively, b) Expected change of nuclear polarization 
(Am), Eq.©. The star indicates the optimal sweep rate /3* 
along the dashed line. 



In the 
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Fourier representation | ip m ) = 
§d9e~ lm6 1 ip$), where \ip m ) and \ip$) are two- 
component spinors V'm) T an d W'Ji V ; e) T '> the 
dynamics for different 9 components decouple. We 
obtain 



(6) 



where H$t + is given in Eq.©. Here the parameter 9 
plays the role of the azimuthal angle of AB nuc t ± in the 
classical problem, Eq.©. Indeed, a state ipg' T oc 5(9— 9 ) 
is a coherent state concentrated near 9q. 

To calculate (Am), we make use of the relation 
m\ tp m ) = 4t § dG 15 (e" lme ) \ipe)- Integrating by parts 
to move the derivative onto | tpe ) , we find 
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where without loss of generality we have taken too = 0. 

To evaluate expression ([7]). we must solve for the two- 
level dynamics of | V'e ) under the time dependent Hamil- 
tonian ([2]). We can write the evolution operator as 



U(9) = 



a 9 b 9 e llpo ^ 
m > e -i Vo {e) a * 



, |a e | 2 + |6 e | 2 = l,(8) 



where <fio(9) = ar g [ve] (see Fig|Tb). With this 
paramctrization, we have ag — a_e, and bg = b-g. 

The initial state | ) = (0, 1) T evolves into \ip g ) — 
(bffe lv °( e \ ag) T . Because ag and bg are even functions of 
9, the only term that contributes to (Am) is that gener- 
ated when the derivative acts on e^W in Eq.([7]). giving 



(Am) 
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where \bg\ 2 is the singlet-triplet transition probability in 
the 9 channel. Comparing this result to that for the net 
electron spin-flip probability, given by P e \ = ^- § \bg\ 2 d9, 



we note that, unlike the situation where only the hyper- 
fine interaction is present, here there is no simple rela- 
tionship between the electron and nuclear spin- flip rates. 

The analysis leading to Eq. © is valid for arbitrary 
time dependence A(t), and in particular can even de- 
scribe cases involving multiple S — T crossings. Remark- 
ably, for very slow sweeps the change in polarization 
becomes sharply quantized. Indeed, since \bg\ 2 — > 1 in 
the adiabatic limit, the integral © is equal to a winding 
number: (Am) = — oz§ d(fo. Thus (Am) = —1 or if 
vg does or does not wind around the origin as 9 traverses 
the Brillouin zone, depending on the relative magnitudes 
of i>hf and uso t 19]. The sign is negative because elec- 
tron transitions from S to T + flip nuclear spins from up 
to down. 

The quantity (Am) exhibits interesting dependence on 
the sweep speed, which can be analyzed most straightfor- 
wardly for linear sweeps A(t) — —fit, when \bg\ 2 is given 
by the Landau-Zener formula ©. Due to the absence of 
dynamics for very fast sweeps, and since (Am) = for 
very slow sweeps in the region \vso\ > |^hf|> the polar- 
ization efficiency is a non-monotonic function of sweep 
rate (see Fig[2jD). In the limit of weak hyperfine interac- 
tion, | «hp | <C I^soIj expanding the exponential in Eq. (J3j , 
wc find 
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(Wt'so)"Ae-\ A = 2Trv 2 so /hf3. (10) 



This expression attains its maximum value |(Am)| max = 
v hf / i ev So) a t the optimal sweep rate /3* = 2ttvq /?i. 

Next we consider a generalization to multiple sweeps 
and show that dynamical polarization is sensitive to nu- 
clear Larmor precession between sweeps. We focus on a 
sequence of two identical gate sweeps, where each individ- 
ual sweep is short on the scale of the nuclear Larmor time, 
but with a waiting time At between sweeps long enough 
to allow nuclear spin precession. Individual sweeps may 
pass through the S — T crossing one or more times. As- 
suming that only one nuclear species contributes to the 
hyperfine field, the change in nuclear polarization due to 
the two sweeps is given by 

/d9 
— (U^dgU), U = U(9 + M)WU (6), (11) 

where U (9) is the evolution matrix for a single sweep, 
given by ©, W is the electronic evolution operator for 
the waiting interval At, the expectation value is taken 
over the state | ) — (0, 1) T , and A9 = lot, At is the 
precession angle, with u>l the nuclear Larmor frequency. 
Due to the periodicity of U(9) in 9, Eq.©, (Am) exhibits 
periodic oscillations in A9 which are a direct manifesta- 
tion of coherent nuclear polarization dynamics. 

The analysis is simple when the detuning |A| is very 
large in the interval between sweeps. Then W ps 
exp(ixcrs), where the phase x is large and very sensi- 
tive to any fluctuations in the waiting time At or other 
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parameters of the system. The resulting randomness 
of x completely suppresses coherence in the electronic 
state between sweeps, while nuclei remain coherent. For 
demonstration, below we choose the particular protocol 
A(t) = f3(t n — t), for \t — t n \ < t, and A = oo otherwise, 
where n = 1, 2, and r <C At = t% — t\ with t\ and t-i the 
times of passage through the avoided crossing. 

For this experimentally relevant case, it is helpful to 
think of the matrix-products defining U and U^ 1 in 
Eq. (fTT|) as sums over amplitudes associated with all pos- 
sible histories of S — T transitions. Fast electron decoher- 
ence between sweeps suppresses the contribution of terms 
arising from non-identical histories in IA and U^ 1 , leaving 
behind a sum of classical probabilities of all trajectories. 
Starting with the singlet initial state and summing over 
all final states, we find 



(Am) = [-(1, 1) d„ § §M(9 + A9)M{9) (0, 1) T ] v 

M(0)J l -^ i 

V ; I p(0)e-*# l-p{9) I 



(12) 



where r\ is an auxiliary variable introduced for book- 
keeping, and p(9) = \bg\ 2 is the transition probability 
in the 9 channel for a single sweep. Here the phase 
<p(6) is the sum of the geometric part <£o(9) — arg[vg], 
the Stokes phase <p s (9), see e.g. Ref. [20| ■ and the dy- 
namical phase <p a d(9) = i E((3t, 9) dt associated with 
adiabatic evolution on a single branch of the two-level 
spectrum of Hamiltonian Hst + , Eq.®, with E(A,9) = 
-^(A/2) 2 + \v B \ 2 . Equation dHJ) then gives: 



(Am) 



where p'{9) = p{9 + A9), ip'(9) = ip{9 + A9). 

The two terms count the contributions from the two 
sweeps. The factor 1 — 2p(9) accounts for the fact that 
a transition from S to T on the first sweep reverses the 
direction of nuclear spin flips in the second sweep. Os- 
cillations and sign reversal of (Am) as a function of the 
precession angle A9 are clearly visible in Fig. |3l where we 
have used direct numerical integration to find the transi- 
tion amplitudes b(9) needed to evaluate expression (p~3]l . 

The dynamical phase contributes to (Am) through 
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where d(p a d/d9 grows logarithmically with the range of 
detunings 2Ao spanned by the sweep. Note that dE/89 
is the group velocity in the Bloch band corresponding to 
the ID quantum walk, Fig. [2ji. 

This contribution describes DNP buildup at times long 
after the level crossing. The in-plane component of elec- 
tron spin maintains a nonzero average value due to the 
presence of the spin-orbit interaction. While this rem- 
nant electron spin polarization can be small, its influence 
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FIG. 3: (Color online) Polarization due to a round trip 
through the S — 7+ crossing with «hf = 0.6, «so = 0.4, Ao = 
50, and precession A8 = cj^At between forward and return 
sweeps, Ea. ()13|) . Here is the nuclear Larmor frequency 
and At is the waiting time. Dashed lines show (Am) calcu- 
lated using ifio{9) and the dynamical phase, Eq. (|14p . neglect- 
ing Stokes phase. The oscillation phase shifts by tt/2 for faster 
sweeps where d ^g d is small (see inset). Units are chosen with 

h=i. 



on nuclear polarization accumulates over a long time, as 
indicated by the 1/(3 dependence. As a result, this con- 
tribution to DNP typically dominates over that of the 
geometric phase <po(9), see Fig El 

Equation (fTTj) can be readily extended to arbitrary se- 
ries of non-identical sweeps. If we take into account clas- 
sical noise in the detuning A(i) during each sweep, the 
primary effect is to flatten the function p{9), causing it 
to saturate toward the value 1/2 for all 9. Because the 
qualitative form of p(9) is preserved, however, oscillations 
resulting from the geometric phase contribution d(fo/d9 
survive under quite general circumstances. 

Resonant effects at the nuclear Larmor frequency also 
can be seen directly in the time dependence of electron 
transport properties. As the transverse polarization pro- 
cesses in time, the matrix element vg traces out the 
dashed circle shown in Fig[lp, leading to a modulation of 
| we | at the nuclear Larmor frequency. This effect can be 
detected as an ac modulation of conductance (current) 
in the system, or as a correlation between the electron 
spin flip probabilities on successive sweeps through the 
avoided crossing. Such measurements would constitute a 
purely electrical detection of nuclear spin dynamics. 

In summary, we have shown that the coherent evolu- 
tion of nuclear spins in quantum dots can be observed 
through oscillations and sign reversal of the nuclear spin 
pumping rate, which occur due to the coherent inter- 
play between hyperfme and spin-orbit couplings. Recent 
observations indicate that this interesting regime is now 
within experimental reach, opening new possibilities to 
explore many-body spin dynamics in the solid state. 
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